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ABSTRACT 



A first-order correction in the a-viscosity parameter of Shakura & Sunyaev has been introduced in the 
standard inviscid and thin accretion disc. A linearised time-dependent perturbative study of the stationary 
solutions of this "quasi-viscous" disc leads to the development of a secular instability on large spatial 
scales. This qualitative feature is equally manifest for two different types of perturbative treatment — a 
standing wave on subsonic scales, as well as a radially propagating wave. Stability of the flow is restored 
when viscosity disappears. 

Subject headings: accretion, accretion disks — hydrodynamics — instabilities — methods: analytical 



1. Introduction 

The compressible, inviscid and thin disc flow has by now become an established model in accretion stud- 
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ies ( Abramowicz & Zurek 1981; Fukue 1987; Chakrabarti 1989; Nakavama & Fukue 1989; Chakrabarti 1990; Kafatos & Yang 
1994 1 Yang & Kafatoslll995||Parievlll996tlMolteni et all 1 9961 ILu et al.lll997tlDaj2002l.lpas et al.ll2003l;]Ravll2003al~ 
Barai et al.ll2004 lDasll2004l kbraham et alfeoodlDas et alj|2006l;lchaudhurv et alhoOftlGoswami et al.ll2007b . This 
is a particularly expedient and simple physical system to study, especially as regards the rotating flow in the innermost 
regions of the disc, in the vicinity of a black hole. Steady global solutions of inviscid axisymmetric accretion on to a 
black hole have been meticulously studied over the years, and at present there exists an extensive body of literature 
devoted to the subject, with especial emphasis on the transonic nature of solutions, the multitransonic character of the 
flow, formation of shocks, and the stability of global solutions under time-dependent linearised perturbations. 

Having stressed the usefulness of the inviscid model among researchers in accretion astrophysics, it must also 
be recognised that this model has its own limitations. It is easy to understand that while the presence of angular 
momentum leads to the formation of an accretion disc in the first place, a physical mechanism must also be found for 
the outward transport of angular momentum, which should then make possible the inward drift of the accreting matter 
into the potential well of the accretor. Viscosity has been known all along to be a just such a physical means to effect 
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infall, although the exact presciption for viscosity in an accretion disc is still a matter of much debate dFrank et al, 
2002). What is well appreciated, however, is that the viscous prescription should be compatible with an enhanc ed 
outward transport of angular momentum. The very well-known a parametrisation of IShakura & Sunvaevl (11973b is 
based on this principle. 

And so it transpires that on global scales — especially on the very largest scales of the disc — the inviscid model 
will encounter difficulties in the face of the fact that without an effective outward transport of angular momentum, 
the accretion process cannot be sustained globally. To address this adverse issue, what is being introduced in this 
paper is the "quasi-viscous" di sc model. This model invol ves prescribing a very small first-order viscous correction 



in the a-viscosity parameter of IShakura & Sunvaevl d 19731) . about the zeroth-order inviscid solution. In doing this, a 



viscous generalisation of the inviscid flow can be logically extended to capture the important physical properties of 
accretion discs on large length scales, without compromising on the fundamentally simple and elegant features of the 



models of axisymmetric flows which involve viscosity (jshakura & Sunvaev 


1973; Liang & Thomson 1980; Pringle 


1981; Matsumoto et al. 1984; Muchotrzeb-Czernv 1986; Abramowicz et al. 


1988; iNaravan & Yilll994 Ichakrabarti 


1996aHdlchen et al.ll997llPeitz & AppJl997tlFrank et albooilAfshordi & Paczvnskil2003lllJmurhan et al.ll2006). 



While many previous works have also taken up the question of the stability of viscous thin d isc accretion dlightman & Eardlev 



1974 IShakura & Sunvaevll 19761: iLivio & Shavivll 19771; lKatolll978t lUmurhan & Shavivll2005l) . the specific objective 
of the present paper is to study the stability of stationary quasi-viscous inflow solutions under the influence of a time- 
dependent and linearised radial perturbation. It has already been shown in some earlier works that inviscid solutions 
are stable under a perturbation of this kind (Ray 2003 at IChaudhury et al.ll2006l) . What has been found through this 



particular work is that with the merest presence of viscosity (i.e. to a first order in a, which itself is much less than 
unity) about the stationary inviscid solutions, instabilities develop exponentially on large length scales. This has dis- 
turbing implications, because all physically meaningful inflow solutions will have to pass through these length scales, 
connecting the outer boundary of the flow with the surface of the accretor (or the event horizon, if the accretor is a 
black hole). 

The perturbative treatment has been executed in two separate ways — as a standing wave and as a high-frequency 
travelling wave, and in both cases the perturbation displays growth behaviour. It need not always be true that standing 
and travelling waves will simultaneously exhibit the same qualitative properties as far as stability is concerned. Many 
instances in fluid d ynamics bear this out. In the case of binary fluids, standing wav es indicate instability as opposed 
to travelling waves ( Cross & Hohenbergll993 : Bhattacharya & Bhattacharjeefl2005|). while the whole physical picture 
is quite the opposite for the fluid dynamical problem of the hydraulic jump dBohr et aL 1993 ; Rav & Bhattachariee 
2004t ISingha et al.l 12005). Contrary to all this the axisymmetric stationary quasi-viscous flow is greatly disturbed 
both by a standing wave and by a travelling wave. This provides convincing evidence of its unstable character, and 
it is very muc h in consonance wit h similar conclusions drawn from some earlier studies. For high-frequency radial 
perturbations lChen & Taaml (1 1993b have found that inertial-acoustic modes are locall y unstable, with a greater degree 
of growth for the outward travelling modes than the inward ones. On the other hand. lKato et al.l (1 1988b have revealed 
a growth in the amplitude of a non-propagating perturbation at the critical point, which, however, stabilises in the 
inviscid regime. 

This kind of instability — one that manifests itself only if some dissipative mech anism (viscous dissip ation in the 
case of the quasi- viscous rotational flow) is operative — is called secular instability ( Chandrasekharfl987 ). It should 
be very much instructive here to furnish a parallel instance of the destabilising influence of v iscous dissipation in a 
system undergoing rotation : that of the effec t of viscous diss i pation in a Maclaurin spheroid dChandrasekharll 1 987b . 
In studying ellipsoidal figures of equilibrium. IChandrasekharl (1 1 987) has discussed that a secular instability develops 
in a Maclaurin spheroid, when the stresses derive from an ordinary viscosity which is defined in terms of a coefficient 
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of kinematic viscosity (as the a parametrisation is for an accretion disc), and when the effects arising from viscous 
dissipation are considered as small perturbations on the inviscid flow, to be taken into account in the first order only. It 
is exactly in this spirit that the "quasi-viscous" approximation has been prescribed for the thin accretion disc, although, 
unlike a Maclaurin spheroid, an astrophysical accretion disc is an open system. 

Curiously enough, the geometry of the fluid flow also seems to be having a bearing on the issue of stability. The 
same kind of study, as has been done here with viscosity i n a rotationa l flow, had also been done earlier for a viscous 
spherically symmetric accreting system. In that treatment dRavll2003bl) viscosity was found to have a stabilising influ- 
ence on the system, causing a viscosity-dependent decay in the amplitude of a linearised standing-wave perturbation. 
This is quite in keeping with the understanding that the respective roles of viscosity are at variance with each other in 
the two distinctly separate cases of spherically symmetric flows and thin disc flows. While viscosity contributes to the 
resistance against infall in the former case, in the latter it aids the infall process. 

Finally, an important aspect of the time-dependent perturbative analysis that may be emphasised is that although 
the flow has been considered to be driven by the Newtonian potential, none of the physical conclusions of this work 
will be qualified in any serious way upon using any of the pseudo-Newtonian potentials ( Paczvnski & Wiita 1980t 



Nowak & Wagoner 199 It lArtemova et alj |1996). which are regularly invoked in accretion-related literature to de- 



scribe rotational flows on to a Schwarzschild black hole, even while preserving the Newtonian construct of space and 
time. This shall be especially true of the flow on large scales, where all pseudo-Newtonian potentials converge to 
the Newtonian limit, and, therefore, the conclusions of this perturbative treatment will also have a similar bearing on 
pseudo-Schwarzschild flows. 



2. The quasi-viscous axisymmetric flow 



For the thin disc, u nder the condition of hydrostatic equilibrium along the vertical direction (IMatsumoto et al 



1984c iFrank et al.ll2.002b . two of the relevant flow variables are the drift velocity, v, and the surface density, E. In 
the thin-disc approximation the latter has been defined by vertically integrating the volume density, p, over the disc 
thickness, H (r). This gives £ = pH, and in terms of S, the continuity equation is set down as 

<9£ 1 d 



dt 



+ -— (£w)=0. 
r or 



(1) 



For a flow driven by the Newtonian potential, V(r) = —GMr^ 1 , assumption of the hydrostatic equilibrium in the 
vertical direction will give the condition H = r(c s /vK), in which the local speed of sound, c s , and the local Keplerian 
velocity, vk, are, respectively, defined as c 2 = jKp" 1 ^ 1 and = GMr^ 1 , with the constants 7 and K deriving 
from the application of a polytropic equation of state, P = Kp 1 , in terms of which, the speed of sound may given as 



Cg = dP/dp. Written explicitly, the disc height is, therefore, expressed as 

1/2 



H = 



GM 



p ( 7 -l)/2 r 3/2 j 



and with the use of this result, the continuity equation could then be recast as 



d 


p (7 + l)/2~ 


1 8 r 


dt 







p (7+l)/2 ur 5/2 



= 0. 



The condition for the balance of specific angular momentum in the flow is given bv (frrank et alj2002b . 

£ (Er 2 fi) +11- UZvr) r 2 Ol = — (^ 
g t \ J ^ r dr LV ) J 2nr \ dr 



(2) 



(3) 



(4) 
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where CI is the local angular velocity of the flow, while the torque is given as 



y = 2iirvY,r \ — — 
or 



(5) 



with v being the kinemat ic viscosity associated with the flow. With the use of the continuity equation, as equation (0 
gives it, and going by the lShakura & Sunvaev ( 1973) prescription for the kinem atic viscosity, v = ac s H, it would be 
easy to reduce equation (0 to the form ( Frank etldl 2002 ; Naravan & Yilll994 ) 



18_ 

v dt 



dr 



2 n) 



d 



pvrH dr 



with fix being defined from vk = rClR. 



Going back to equation (0, a new variable is defined as / 



apHc^r 



(7+l)/2, 



dfl 
dr 



(6) 



- D / 2 , whose steady value, as it is very 



easy to see from equation 0, can be closely identified with the constant matter flux rate. In terms of this new variable, 
equation can be modified as 



d_ 

dt 



(7+l)/2 



J_d£ 

r-5/2 Q r 



0, 



while equation (O can be rendered as 

v ot ' 



/ jK \ ld_ 
a \GMj l~dr 



p2y3 



dr 



(7) 



(8) 



The inviscid disc model is given by the requirement that r 2 £! = A, in which A is the constant specific angular mo- 
mentu m. The quasi-viscous di sc that is being proposed here will introduce a first-order correction in terms involving a, 
the lShakura & Sunyaevl (119731) viscosity parameter, about the constant angular momentum solution. Mathematically 
this will be represented by the prescription of an effective specific angular momentum, 



A eff (r) = r 2 Q = A + ar 2 tt, 



(9) 



with the form of £1 having to be determined from equation ©, under the stipulation that the dimensionless a-viscosity 
parameter is much smaller than unity. This smallness of the quasi-viscous correction induces only very small changes 
on the constant angular momentum background, and, therefore, neglecting all orders of a higher than the first, and 
ignoring any exp licit time-variation of the vi s cous correction term, the latter being a standard method ado pted also for 
Keplerian flows jLightman & Eardleyl[l974l Ishakura & Sunyaevl[l976[ |Pringldfl98lt iFrank et al.ll2002l) . the depen- 
dence of fi on v and p is obtained as 



n = 



2A 



GM 



/ 2 »k 

p2 v 3 r 3 



/ 2 ^K 
p2 v 3 r 3 



ldf 

f dr 



dr 



(10) 



The effect of the absence and the presence of a small viscous correction to the inviscid background flow has been 
schematically shown in Figs.Q]&[2l respectively. 



The stationary solution of equation (0 can be easily obtained as a first integral, and with the help of equation (0 
it can be set down as 



2n 



GM 



1/2 



p (7+l)/2 ur 5/2 = 



(ID 



with m being the conserved matter inflow rate. The negative sign arises because for inflows, v goes with a negative 
sign. Further, under stationary conditions, equation (0 can be written in a modified form as 



Aeff(r) = A - 2a A 



(A. 

\vv K 



(12) 



Fig. 1 . — A schematic representation of the axisymmetric inviscid disc. Viewed from the top (along the vertical axis 
of the disc) the flow will be seen to form closed circular paths. 




Fig. 2. — A schematic representation to trace the spiralling behaviour of the physical flow trajectories upon the 
introduction of a small viscous correction to the background inviscid flow. 
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From equation ( fTTT i. with p approaching a constant asymptotic value on large length scales, the drift velocity, v, can 
be seen to go asymptotically as r~ 5 / 2 . Bearing in mind that for inflows, v < 0, the asymptotic dependence of the 
effective angular momentum can be shown to be 



A e fi(r) ~ A + 2a A 



(13) 



with r s being a scale of length, which, to an order-of-magnitude, is given by r 3 ~ GMm[c^ (oo)p(oo)] -1 . This 
asymptotic behaviour is entirely to be expected, because the physical role of viscosity is to transport angular momen- 
tum to large length scales of the accretion disc. 

Lastly, the equation for radial momentum balance in the flow will also have to be modified under the condition 
of quasi-viscous dissipation. This has to be done according to the scheme outlined in equation (O by which, the 
centrifugal term, A 2 ff (r)/r 3 , of the radial momentum balance equation, will have to be corrected upto a first order in 
a. This will finally lead to the result 



dv dv IdP „„ . A 2 A / , 
^7 +v— + -— + V'(r) - — - 2a^ ( r 2 tt 
at or p or r 6 r A 



0, 



(14) 



with Cl being given by equation ([Tol l, and P being expressed as a function of p with the help of a polytropic equation 
of state, as has been mentioned earlier. The steady solution of equation ( fT4b is given as 



dv 1 dP T . . , 
v— + -— : + V'(r) 
ar p ar 



- — +4a— -2_ =0, 

r A r^ \ tok / 



(15) 



whose first integral cannot be obtained analytically because of the a-dependent term. In the inviscid limit, though, 
the integral is easily obtained. T his case will be governed by conserved cond itions, and its solutions have been 



well-known in accretion literature dchakrabartill 1 98% IPashooi 



Dasetal 



20031) . They will either be open solutions 



passing through saddle points or closed paths about centre-type points. The slightest presence of viscous dissipation, 
however, will radically alter the nature of solutions seen in the inviscid limit, and it may be easily understood that 
solutions formi ng closed paths about centre-type points, will, under conditions of small- visco us correction, be of the 
spiralling kind ( Liang & Thomson! 1980t Matsumoto et aljfl984 ; Afshordi & Paczvnskill2003 ). This state of affairs is 
appreciated very easily by the analogy of the simple harmonic oscillator. In the undamped state the phase trajectories 
of the oscillator will, very much like the solutions of the inviscid flow, be either closed paths about centre-type points 
or open paths through saddle points. With the presence of even very weak damping the closed paths change into 
spiralling solutions. 

After having understood the qualitative nature of the stationary flows in the quasi-viscous accretion disc, as given 
by equations ( fTTT l, ( fT2] i and ( fT5l ), it will now become possible to carry out a real-time linear stability analysis about the 
stationary solutions, with the help of equations (O and ( TBI ). 



3. An equation for a time-dependent perturbation on stationary solutions 

About the stationary solutions of the flow variables, v and p, a time-dependent perturbation is introduced accord- 
ing to the scheme, v(r, t) = vo(r) + v'(r, t), p(r, t) = po(r) + p'(r, t) and /(r, t) = fo(r) + f'(r, t), in all of which, 
the subscript "0" implies stationary values, with fo especially, as can be seen from equation ([3]), being a constant. 
This constant, as it is immediately evident from a look at equation ( fTTT l. is very much connected to the matter flow 
rate, and, therefore, the perturbation /' is to be seen as a disturbance on the steady, constant background accretion 
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rate . For spherically symm etric flows, this Eulerian perturbation scheme has been applied by iPetterson et al. I dl98d) 
and Theuns & David ( 1992), while for invisc id axisymmetric flows, the same method has been used equally effectively 
bv lRavl d2003albl) and lchaudhurv et aD d2006l) . 



The definition of / will lead to a linearised dependence among /', v' and p' as 



L 

fo 



7 + 1 



P_ 

Po 



V 

VQ ' 



while from equation ©, an exclusive dependence of p' on /' will be obtained as 



df/_ 
~dt 



■0 



2 V 0P0 

fo 



d_r_ 

dr 



with (3 2 = 2(7 + 1) 1 . Combining equations ( fT6l ) and ( fT71 ) will render the velocity fluctuations as 



dt 



vo 
fo 



df_ 

dt 



Vo- 



d_l_ 

dr 



which, upon a further partial differentiation with respect to time, will give 



i:) 2 



dt 2 



vo 
fo 



dt 2 



d_ 

dr 



df_ 

dt 



From equation dT4b the linearised fluctuating part could be extracted as 



d_ 

dr 



v v 



Po 



AaX 1 



2^-2^-3^ 



d_ 

dr 



dr 



0. 



(16) 



(17) 



(18) 



(19) 



(20) 



fo Po v 

in which a = c 2 a / (vqvk) and c s o is the local speed of sound in the steady state. Differentiating equation (|20| > partially 
with respect to t, and making use of equations (1171 . ( fT8l and (fT9l to substitute for all the first and second-order 
derivatives of v' and p', will deliver the result 



d^r_ 

dt 2 



dr 





1 


d ' 




«0 


dr 



VO [Vq 



P 2 clo 



df_ 

dr 



AaX 1 



vor° 



d_r_ 

dt 



/3 7 -l\ df 



a 



V 7+ 1 

d 



vo 



dr 



dr 



d_r_ 

dt 



dr 



= (21) 



entirely in terms of /'. This is the equation of motion for a perturbation imposed on the constant mass flux rate, fo, 
and it shall be important to note here that the choice of a driving potential, Newtonian or pseudo-Newtonian, has no 
explicit bearing on the form of the equation. 



4. Linear stability analysis of stationary solutions 

With a linearised equation of motion for the perturbation having been derived, a solution of the form /' (r, t) = 
gui(r) exp(— iut) is applied to it. From equation (f21~b . this will give 

2 „• d , n Id 

u g u + 2iw— (voffcj v- 

ar Vq dr 

+7 /-(£)*]-•■ <22 » 

The perturbation may now be treated as a standing wave spatially confined between two suitably chosen boundaries, 
as well as a radially propagating high-frequency wave. These two distinct cases will be taken up separately in what 
follows, to see how stationary solutions are affected by the perturbation. 



Vo («o 



B 2 2 \ ^ 



4aA 2 



Vor a 



( 37 - 1 \ dg^ 
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4.1. Standing waves 



It can be easily appreciated that with vis cous dissipation present in the accreting system, multiply- valued solutions 
in the phase portrait are a distinct possibility ( Liang & Thomson 1980t Abramowicz & Katoll989 ; Afshordi & Paczvnski 
2003). This, however, is physically not feasible for a fluid flow, and, therefore, for all solutions which are multiply- 
valued about a critical point, it should be necessary to have the inner bran ch of a solution f itted to its outer branch via 
a shock, which is a standard practice in general fluid dynamical studies (Bo hr et alj|l993l) . The outer branch of this 
discontinuous solution will connect the shock with the outer boundary of the disc itself, and in the phase portrait of the 
flow, there will exist an entire family of these outer solutions which arguably shall be subsonic. It would be worthwhile 
at this point to remember that for this kind of a thin disc system, the local speed of acoustic propagation would go as 
Peso, and subsonic solutions would have their bulk flow velocity, vq, at less than this value. Over the entire subsonic 
range, at two chosen points, the perturbation may be spatially constrained by requiring it to be a standing wave, which 
dies out at the two chosen boundaries. The outer point could suitably be chosen to be at the outer boundary of the flow 
itself, where, by virtue of the boundary condition on the steady flow, the perturbation would naturally decay out. The 
inner boundary of the standing wave is to be chosen infinitesimally close to the shock front, through which the flow 
will be discontinuous, and the perturbation will be made to die out in its neighbourhood. Between these two points, for 
solutions which are entirely subsonic, it should be necessary to multiply equation (l22l by vog^, and then carry out an 
integration by parts. The requirement that all integrated "surface terms" vanish at the two boundaries of the standing 
wave, will give a quadratic dispersion relation in u), which will read as 



Auj 2 +Buj + C = Q, 



(23) 



in which the coefficients of each successive term will be given by 

A = J (v gl) dr, B = -4k*A 2 

and 



dr 



dr + 2a\ 2 (37 - 1) 



dr 



dr. 



It is possible now to find a solution for w, but of greater immediate significance is the fact that this solution will 
yield an a-dependent real part of the temporal component of the perturbation, which will go as 



3? (-iw) = 2a 



with f(r) itself being expressed as 

S(r) = 



A 2 



vog^r 



(v gl) dr 



(24) 



(25) 



in which the Mach number, M. = vq/c s q. It is globally true that «k ~ r^ 1 / 2 . In this situation, under the assumption 
that the Mach number has a power-law dependence on the radial distance, given as M 2 ~ r~ e , it is possible to write 
down, 



dQn|) _ 5 
d(lnr) ~ e ~ 2' 



(26) 



On large length scales, with c s q approaching a constant ambient value, and with vq ~ r 5 / 2 from the continuity 
condition, it is easily deduced that £(r) ~ r 5 / 2 . This indicates that on large length scales, the amplitude of the spatially 
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constrained standing-wave perturbation grows in time, i.e. the subsonic solutions on which this perturbation has been 
imposed, display unstable behaviour. It is easy to see that when a vanishes, all stationary solutions restricted by the 
condition vo < (3c s q, are stable under a standing-wave perturbation, which displays a pure l y oscil latory behaviour 
with no growth in amplitude. For the inviscid di sc (with a = 0) th is has been shown by Ray (2003a). Regarding this 
point, it will be instructive to mention again that lKato et al.l d 1988b have pointed to the existence of a non-propagating 
growing perturbation localised at the critical point, but interestingly enough they have also found that this instability 
disappears in inviscid transonic flows. 



4.2. Radially propagating waves 

The perturbation is now made to behave in the manner of a radially travelling wave, whose wavelength is suitably 
constrained to be small, i.e. it is to be smaller than any characteristic length scale i n the system. A pert urbative 
treatment of this na ture has been carried out before on spherically symmetric flows dPetterson et al. 198dh and on 
axisymmetric flows (Ray 2003a ; Chaudhurv et al.lBo06l) . In both these cases the radius of the accretor was chosen as 
the characteristic length scale in question, and the wavelength of the perturbation was required to be much smaller than 
this length scale. In this study of an axisymmetric flow driven by a Newtonian potential, the radius of the accreting 
star, r*, could be a choice for such a length scale. As a result, the frequency, cu, of the waves should be large. 

An algebraic rearrangement of terms in equation ( 1221 will lead to an integro-differential equation of the form 

r^TT + 2^ - + T / a- (t) dr = °> (27) 

dr z dr J \dr J 

with its coefficients being given by 

V = vl - /? 2 c s 2 , Q = 3v Q ^ - 1 A ( Uo /3 2 c 2 ) - 2icuv - 2aA 2 (3 7 - 1) ^f, 

TL = 2iuj—^-+lu 2 and T= — . 

dr vqt 6 

At this stage, bearing in mind the constraint that lu is large, the spatial part of the perturbation, (r), is prescribed 
a s 9uj(r) = exp(s), with the function s itself being represented as a power series of the form 



<r) V (28) 



k n (r) 



n=-l 

The integral term in equation d27] i. can, through some suitable algebraic substitutions, be recast as 

J 3<o (^) dr = J exp(s) ds = g„{r)S, 

with S itself being given by another power series as 

^Enr 1 ^. (29) 

m— 1 

Following this, all the terms in equation ( |27] | can be expanded with the help of the power series for g^ir). Under the 
assumption (whose self-consistency will be justified soon) that to a leading order 

d<7 da f dfe_i \ 1 
S ~ — — — — — I cu- 



ds dr V dr 
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the three successive highest-order terms (in a decreasing order) involving u> will be obtained as uj 2 , uj and u . The 
coefficients of each of these terms are to be collected first and then individually summed up. This is to be followed by 
setting each of these sums separately to zero, which will yield for uj 2 , uj and uj°, respectively, the conditions 



(«8 



P 2 cf 



dk-i 
dr 



d 2 fc_i dfc_i dfc 



dr 2 



dr dr 



3v 



dr 



- 2iu 

1 d 

Vq dr 



dfc_i 
dr 



1 = 



(30) 



{v (3 2 c 2 ) - 2a\ 2 (37 - 1) 



(3 2 a 



dfc_ 



dr 



dfc dv 
2ivq— : 2i— — = (31) 



dr 



dr 



and 



d 2 fc dfc_i dfci 



dr 2 



dr dr 









+ 



3v Q 



dr 



1 d 

vq dr 



(v (3 2 c 2 )~2aX 2 (37-I) 



/? 2 a 



dfc 
dr 



2i» 



dfci 



4iaA 2 



1 dcr / dfc_i 



dr vor 3 dr \ dr 

Out of these, the first two, i.e. equations (f30b and OTb . will deliver the solutions 



V ± /5c s0 



dr 



0. (32) 



(33) 



and 



respectively. 



fc = - - In {v (3c s0 ) ± a\ 2 (3j - 



(3c s0 {vq ± /3c s0 ) 
u u K r 3 (w 2 - /3 2 c 2 



■ dr, 



(34) 



The two foregoing expressions give the leading terms in the power series of <? w (r). While dwelling on this 
matter, it will also be necessary to show that all successive terms of s(r) will self-consistently follow the condition 
w _n |A; n (r)| ^> cj _ (" +1 ^|A;„ + i(r)|, i.e. the power series given by g u (r) will converge very quickly with increasing 
n. In the inviscid limit, this requirement can be shown to be very much true, considering the behaviour of the first 
three terms in k„ (r) from equations d33l . (1341 and d32l . These terms can be shown to go asymptotically as k-\ ~ r, 
k ~ lnr and fej ~ r _1 , given the condition that w ~ r -5 / 2 on large length scales, while c s q approaches its 
constant ambient value. With the inclusion of viscosity as a physical effect, it can be seen from equations (|33l 
and (l34l . respectively, that while remains unaffected, ko acquires an a-dependent term that goes asymptotically 
as r. This in itself is an indication of the extent to which viscosity might alter the inviscid conditions. However, 
since a has been chosen to be very much less than unity, and since the wavelength of the travelling waves is also very 
small, implying that w ^> (vq ± /3c s o)/r^, the self-consistency requirement still holds. Therefore, as far as gaining 
a qualitative understanding of the effect of viscosity is concerned, it should be quite sufficient to truncate the power 
series expansion of s(r) after considering the two leading terms only, and with the help of these two, an expression for 
the perturbation may be set down as 



f'(r,t) 



exp 



±a\ 2 (37 - I] 



Pc s0 (v ± /3c s0 ) 
v v K r 3 {vl - f3 2 c 2 ) 



dr 



exp 



V ± /3c s0 



■ dr 



(35) 



which should be seen as a linear superposition of two waves with arbitrary constants A + and A- . Both of these two 
waves move with a velocity (3c s q relative to the fluid, one against the bulk flow and the other along with it, while 
the bulk flow itself has a velocity vq. It should be immediately evident that all questions pertaining to the growth 
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or decay in the amplitude of the perturbation will be crucially decided by the real terms delivered from fco. The 
viscosity-dependent term is especially crucial in this regard. For the choice of the lower sign in the real part of /' in 
equation d35l >, i.e. for the outgoing mode of the travelling wave solution, it can be seen that the presence of viscosity 
causes the amplitude of the perturbation to diverge exponentially on large length scales, where c s o — c s (oo) and 
vo ~ r -5 / 2 , with ~vq being positive for inflows. The inwardly travelling mode also displays similar behaviour, albeit 
to a quantitatively lesser degree. It is an easy exercise to see that stability in the system would be re stored for the limit 
of a = 0, and this particular issue has been discussed by iRayl d2003al) and Chaudhury et al.1 (hood) . The exponential 



growth behaviour of the amplitude of the pertu rbation, therefore, is exclusively linked to the presence of viscosity. 
Going back to a work of IChen & Taaml (11993b . it can be seen that the inertial-acoustic modes of short wavelength 
radial perturbations are locally unstable throughout the disc, with the outward travelling modes growing faster than 
the inward travelling modes in most regions of the disc, all of which is very much in keeping with what equation ( |35t 
indicates here. 



With the help of equation (fTTT i it should be easy to express the density fluctuations in terms of /' as 



p \iudrj f '" 

and likewise, the velocity fluctuations may be set down from equation ( [ToT l as 

v \ iujdr) f ' 



(36) 



(37) 



In a unit volume of the fluid, the kinetic energy content is 



1 2 
£kin = - Oo + P ) Of) + V ) , 



(38) 



while the potential energy per unit volume of the fluid is the sum of the gravitational energy, the rotational energy and 
the internal energy. For a quasi-viscous disc, to a first order in a, this sum is given by 



pot 



V(r) 



\ 2 

A cff 

2r 2 



d 
dp 



r a 2 

dpi 



2 (Poe) 



(39) 



where e is the internal energy per unit mass ( Landau & Lifshit3ll987 ). In equation d39l the effective angular momen- 
tum for the quasi-viscous disc will have to be set up as a first-order correction about the inviscid conditions. Following 
this, a time-dependent perturbation has to be imposed about the stationary values of v and p. All first-order terms 
involving time-dependence in equations d38l ) and ( [391 will vanish on time-averaging. In this situation the leading 
contribution to the total energy in the perturbation comes from the second-order terms, which are all summed as 



1 /2 . , , . 1 ,2 

pert = -^POV +V p V + -p 



d 2 

dpi 

1 



2 (A>e) 



2aA 2 



-6 



-2 



PT_ 

Poh 



3 



v'p' 
vopo 



-6 



f'v' 

foVQ 



dr 



(40) 



,Po) J a dr (/o) a \° (/ 2 p ^ Vq) dr (/ 0/ 

In the preceding expression all terms involving p' and v' can be written in terms of /' with the help of equations < 
and d37l i. in both of which, to a leading order, s ~ ojk-x. This is to be followed by a time-averaging over f' 2 , 
which will contribute a factor of 1/2. The total energy flux in the perturbation is obtained by multiplying £ pex t by the 
propagation velocity (vq ± (3c s q) and then by integrating over the area of the cylindrical face of the accretion disc, 
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which is 2itrH. Under the thin-disc approximation, H <C r, this will make it possible to derive an estimate for the 
energy flux as 

l-/3 2 (2-/,)- 



Hr) 



irP 2 A 2 ± 



fo 



GM 



±1 



2/3 (M± (3) 



1 - 



2a\ 2 ^p 
(3 2 M 2 v v K r 2 



x exp 

in which M. is the Mach number, as defined earlier, while 



±2a\ 2 (3 7 -l) J 



(3c s0 (vp ± /3c s0 ) 
v v K r 3 (v$ - 2 c 2 Q ) 



dr 



(41) 



po 



1 sO 



d 2 (ppe) 
dpi 



and 



(p 2 -1) 2 M 2 ±M(3 (/3 2 - 4) + (i 2 [1 ± 2/3M + (3 2 p 



When Ai — ► on large length scales, ip converges to a finite value. However, on these same length scales, what will 
not converge are the two terms involving a in equation (PTTT i. Under the asymptotic conditions on vq and c s q, discussed 
earlier, one term will diverge exponentially as r, while another will have a power-law growth behaviour of r 6 . The 
quasi-viscous disc will, therefore, be most pronouncedly unstable on large scales under the passage of a linearised 
radially propagating high-frequency perturbation. It is easy to check that under inviscid con ditions, with a = 0, and 
for an adiabatic perturbation with p = l, the disc will immediately revert to stable behaviour jRayll2003al) . 



5. Concluding remarks 

The instability of the quasi-viscous disc raises some doubts, the primary one of these being on the possibility 
of a long-time evolution of the disc towards a stationary end. The quasi-viscous disc being a dissipative system, 
i.e. its energy being allowed to be drained away from this system, there cannot be any occasion to look for the 
selection of a particular solution, and a selection criterion thereof, on the basis of energy minimisation, as it can 



be done for an idealised in viscid flow, such as the Bondi solution in spherical symmetry ( Bondi 1952 ; Garlick 19791; 



Rav & Bhattacharieel2002l) . This, of course, shall also affect the flow rate, which, in this treatment, has been perturbed 



and has been found unstable. As a result of all this the whole system has been left without any well-defined criterion 
by which it could guide itself towards a steady state, transonic or otherwise. 

One very important physical role of viscosity in an accretion disc is that it determines the distribution of matter in 
the disc. The manner in which viscosity redistributes an annulus of matter in a Keplerian flow around an accretor is very 
well known, with the inner region of this disc system drifting in because of dissipation, and consequently, through the 
conservation of angular momentum and its o utward transport, making it nece ssary for the outer regions of the matter 



distribution to spread even further outwards dPringlelll98ll Fran k et al ] l2002b . This state of affairs is qualitatively not 



altered in anyway for the quasi-viscous flow, except for the fact that with viscosity being very weak in this case, the 
outward transport of angular momentum can be conspicuous only on very large scales. It may rightly be conjectured 
that the instability that develops on the large subsonic scales of a quasi-viscous disc is intimately connected with the 
cumulative transfer of angular momentum on these very length scales. The accumulation of angular momentum in 
this region may create an abrupt centrifugal barrier against any further smooth inflow of matter. This adverse effect, 
on the other hand, could disappear if there should be some other means of transporting angular momentum from the 
inner regions of the disc. Astrophysical jets could readily afford such a means, insofar as jets actually cause a physical 
drift of angular momentum vertically away from the plane of the disc, instead of along it. 
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